Junctions of one-dimensional quantum wires — correlation effects in transport 
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We investigate transport of spinless fermions through a single site dot junction of M one- 
dimensional quantum wires. The semi-infinite wires are described by a tight-binding model. Each 
wire consists of two parts: the non-interacting leads and a region of finite extent in which the 
fermions interact via a nearest-neighbor interaction. The functional renormalization group method 
is used to determine the flow of the linear conductance as a function of a low-energy cutoff for a 
wide range of parameters. Several fixed points are identified and their stability is analyzed. We 
determine the scaling exponents governing the low-energy physics close to the fixed points. Some 
of our results can already be derived using the non-self-consistent Hartree-Fock approximation. 



I. INTRODUCTION 



In one spatial dimension correlation effects strongly 
influence the low-energy physics of many-fermion sys- 
tems. Such systems cannot be described as Fermi liquids, 
but are classified as Tomonaga-Luttinger liquids (TLLs), 
which are characterized by a vanishing quasi-particle 
weight and power-law scaling of correlation functions. 1 
For spin rotationally invariant interactions and spinless 
fermions, on which we focus here, the exponents of the 
different correlation functions can be expressed in terms 
of a single number, the TLL parameter K . It depends 
on the parameters of the chosen model, in particular the 
strength of the two-particle interaction. For vanishing 
interaction K = 1, while < K < 1 for repulsive inter- 
action and K > 1 in the attractive case. As indicated by 
the singular behavior of the density response function at 
momentum transfer 2kp, 2 ' 3 with kp the Fermi momen- 
tum, a TLL reacts quite differently to an inhomogeneity 
than a Fermi liquid. The physics of inhomogeneous TLLs 
can conveniently be studied investigating transport prop- 
erties. 

The simplest junction is a single impurity in an infi- 
nite TLL wire. The transport through such a system has 
intensively been studied in the past. Using the renormal- 
ization group (RG) language the single impurity problem 
can be characterized by two fixed points (FPs). 4 One 
is the "perfect chain" FP at which the impurity effec- 
tively vanishes and the conductance takes its maximal 
value. For TLL wires that are "smoothly" coupled to 
non-interacting leads, a situation we consider here, the 
latter is given by G = e 2 /h. 5 The correction to the FP 
conductance asymptotically scales as s 2 ^ -1 ^, with s be- 
ing the largest (but still asymptotically small) energy 
scale (e.g. temperature, bias voltage, external infrared 
cutoff). For < K < 1 the exponent is negative and the 
FP is unstable, while it is stable for K > 1. The other 
FP is the "decoupled chain" FP at which G = 0. The 
correction scales as s 2aB , with as — 1/K — 1. The FP is 
stable for repulsive interactions and unstable in the at- 
tractive case. The exponent as characterizes the power- 
law behavior of the local one-particle spectral weight of 
a TLL with an open boundary close to the boundary. 4,6 



The flow from one to the other FP is described by a in- 
dependent one-parameter scaling function. The scaling 
behavior of the conductance has been demonstrated for 
a simplified effective low-energy model 4,7-10 as well as for 
a microscopic lattice model. 11,12 

Recently single-walled carbon nanotubes were used 
to experimentally realize junctions of several quasi one- 
dimensional quantum wires. 13,14 They might form the 
basis of future nano-electronic devices. Taking into ac- 
count the fermion interaction, models for different types 
of junctions and networks of TLLs have been investigated 
theoretically using a variety of methods. 15-22 These stud- 
ies left open several interesting questions. Already the 
low-energy physics of the three wire Y-junction is much 
richer than that of the single impurity problem. 

We here study the transport through a single site dot 
junction of M = 2, 3, . . . semi-infinite wires, each de- 
scribed by a microscopic lattice model, at temperature 
T = 0. To obtain the conductance between the M legs we 
mainly use an approximate technique that is based on the 
functional renormalization group (fRG) method. 23-25 It 
has earlier been successfully applied to describe the trans- 
port in a TLL with a single impurity 11,12 and a double 
barrier, 12,26 the latter allowing for resonant tunneling. 
The approximations lead to reliable results for not too 
strong interactions with TLL parameter 1/2 < K < 3/2. 
In particular, for a single impurity the power-law scaling 
of the conductance discussed above is reproduced with 
exponents that agree with the exact ones to leading order 
in the interaction. For the M-leg junction we investigate 
the RG flow for a wide range of parameters and identify 
the FPs. We numerically determine the exponents of the 
power-law corrections to the FP conductances that gov- 
ern the low-energy physics close to the FPs. They depend 
on the interaction and the number of wires M. Most of 
these exponents have not been determined before. As 
in our approximation terms of second order in the inter- 
action are only partly included the exponents can only 
expected to be correct to leading order in the interac- 
tion. In a short publication we have earlier verified that 
for a specific type of triangular three wire junction (not 
discussed in the present publication), for which an ex- 
act result is available, 19 we indeed reproduce the scaling 
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exponent to leading order. 22 

For a specific set of junction parameters the fRG study 
is supplemented by results for the conductance obtained 
using the non-self-consistent Hartree-Fock approxima- 
tion (HFA) and Fermi's golden rule like arguments. The 
HFA allows us to analytically calculate one of the scaling 
exponents. It has earlier been shown that this approxi- 
mation leads to meaningful results for the power-law scal- 
ing of the one-particle spectral weight in a TLL with an 
open boundary. 27 

The paper is organized as follows. In Sect. II we intro- 
duce our model of the M-wire junction. The fRG based 
approximation scheme is discussed in Sect. III. Using 
single-particle scattering theory in this section we also 
derive equations relating the conductance to matrix ele- 
ments of an auxiliary Green function and the dot Green 
function. They can be used to reduce the numerical effort 
for solving the fRG flow equations and to gain a deeper 
understanding of our findings for the conductance. In 
Sect. IV we apply the HFA to determine scaling expo- 
nents for a certain class of symmetric junctions. Our 
fRG results for the FPs, the scaling exponents of the cor- 
rections to the FP conductances, and the general RG flow 
are presented in Sect. V. We conclude with a summary 
and an outlook in Sect. VI. 



II. THE MODEL 

Each of the M quantum wires that meet at that sin- 
gle site dot junction is described by the lattice model 
of spinless fermions with nearest-neighbor hopping. The 
semi-infinite wires can be divided in two sections: the 
lead with lattice sites j > N in which the fermions are 
assumed to be non-interacting and the interacting wire 
with nearest-neighbor interaction across the bonds of the 
sites j G [1, N]. Fig. 1 shows a sketch of our system. We 
here focus on the half-filled band case. The results are 
generic also for other fillings. 

The Hamiltonian reads 



H — -ffkin + flint + -Hjunc ■ 

The kinetic energy is modeled by 



(1) 



M oo 

ffkin = ~* X X ( C }+1," C J> + C L C 3+h" ) ( 2 ) 
u=l j = l 

where we used standard second-quantized notation with 
Cj v and c- v being creation and annihilation operators on 
site j of wire v, respectively. From now on we set t = 1, 
i.e. measure energies in units of t. 

As the part of the Hamiltonian containing the interac- 
tion we take 

M N-l 

H ^ = X X Ujj+i - V2] [n j+ i, v - 1/2] , (3) 




FIG. 1: A single site dot junction of M quantum wires. Across 
the bonds of the lattice sites j = 1, . . . , N (small filled circles) 
the fermions interact via a nearest-neighbor interaction, while 
they are non- interacting in the leads with j > N (solid line). 
The hopping amplitude between the first site of wire v = 
1, . . . , M and the dot site is t v . The on-site energy on the dot 
site (large filled circle) is V . 



with the local density n j v 



The interaction 



Ujj+i is assumed to be independent of the wire index and 
acts only between the bonds of the sites 1 to V, that de- 
fine the interacting wire. Within this region it is allowed 
to depend on the position. By subtracting the average 
filling 1/2 from the density n^ v we prevent a depletion 
of the interacting part of the wire. The chemical poten- 
tial corresponding to half- filling is /x = 0. To avoid any 
fermion backscattering at the contact between the lead 
and the interacting wire, Ujj+i is turned on smoothly 5 
starting at zero across the bond (N, N + l) and approach- 
ing its bulk value U at bond (N — j s , N — j s + l). 11 . 12 ,28 
More explicitly we use 



U, 



U i<> 



arctan [s n(2[j — N] + j s )/js] 
arctan[s tt] 



(4) 



for j = N - j s , ...,N and U jij+1 = U for 1 < j < N - 
j s . The larger N the smoother the interaction has to be 
switched on. We here consider interacting wires of up to 
N = 10 5 sites for which j s = 32 and s = 2 turned out to 
be sufficient. For these parameters the backscattering at 
the lead-interacting wire contact is less than 10~ 4 % and 
can thus be neglected. The results do not dependent on 
the detailed shape of the envelope function as long as it 
is sufficiently smooth. 

The model corresponding to the Hamiltonian i/kin + 
H[ n t with interaction U across all bonds (not only the 
ones within [1, Nj) and M = 1 shows TLL behavior for 
\U\ < 2 with a TLL parameter (for half-filling) 29 



K = 



2 / U 
— arccos 

TT V 2 



-l -1 



To leading order in the interaction it is given by 
K = 1-- + 0{U 2 ) , 

7T 

an expression we repeatedly refer to further down. 



(5) 



(6) 
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The junction we model by 



M 



Y,t v {c{ v d +Sc l<v )+VSd , (7) 



with d) and d being creation and annihilation opera- 
tors on the dot site, respectively. It is parameterized by 
the hopping amplitudes t v > connecting the wire v to 
the dot and the on-site energy V > on the dot. For 
M = 2 the junction is equivalent to a local impurity in 
an infinite wire. Applying the fRG for this case we re- 
cover the results for the conductance obtained earlier (see 
below). 11 - 12 

Note that in our Hamilltonian the fermions on the dot 
site do not interact with the fermions on the first lattice 
sites of the wires. Including such additional interactions 
does not lead to any changes of the FP structure and 
scaling exponents investigated here, as we have verified 
explicitly. We exclude such terms from our model as 
otherwise we later would have to introduce renormalized 
junction parameters which would lead to an unnecessary 
proliferation of symbols. 12 



III. THE METHOD 

At T = all inelastic processes are frozen out and the 
linear conductance G UM i from wire v to wire v' can ex- 
actly be expressed in terms of a real space matrix element 
of the one-particle Green function Q(e + zO) evaluated at 

e = n = O, 12 - 30 



\t vy \ 2 =l\{N,v\G{Q + M)\N, v' 



(8) 



Here \N, v) denotes the Wannier state centered on site 
TV of wire v and \t v y | is the effective transmission from 
wire v to is', with v ^ v 1 and v,v' E [1,M]. Note that 
Q(e + iO) must be calculated in the presence of the non- 
interacting leads, the junction, and the interaction. 



A. The functional renormalization group 

To obtain an approximation for the Green function we 
use the fRG. A detailed account of the method was given 
in Refs. 12 and 31. We here only present the approximate 
flow equations (which are then integrated numerically), 
describe the most important steps to derive them, and 
give details specific to the present junction geometry. 

An infrared cutoff A is introduced by replacing the 
non-interacting imaginary frequency propagator Q of the 
system by the A dependent propagator 



0oM = 6(M-A)0oM 



(9) 



The cutoff runs from A = oo down to A = 0, at which 
Q A=0 (iuj) = Q (iu) is reached and the cutoff- free prob- 
lem is recovered. Using the generating functional for 



one-particle irreducible vertex functions, with Qq as the 
non-interacting propagator, an infinite hierarchy of cou- 
pled flow equations for the self-energy, the effective two- 
particle interaction, and higher order vertex functions 
is derived. 23-25 It is truncated by neglecting the three- 
particle vertex, which is a valid approximation as long as 
the two-particle vertex does not become too large. The 
two-particle vertex projected on the Fermi points is pa- 
rameterized by an effective nearest-neighbor interaction 
U A . The flow equation for the latter is obtained by con- 
sidering a single infinite chain with interaction across all 
bonds and neglecting self-energy corrections. 31 It can be 
integrated and at half-filling U A is given by 



u 



A - U/(2n) 



1 



(10) 



The A-dependent two-particle vertex is then approxi- 
mated by a frequency independent nearest-neighbor in- 
teraction of strength U A . In the case where the inter- 
action depends on position, as an additional approxima- 
tion we apply Eq. (10) locally for each bond. As a con- 
sequence of the assumed frequency independence of the 
effective interaction also the self-energy does not depend 
on lu. In the exact solution an ^-dependence is generated 
to order U 2 (bulk TLL behavior). This exemplifies that 
in our approximation for the self-energy terms of order 
U 2 are only partly included. 

With these approximations the self-energy is diagonal 
in the wire index v and tridiagonal in the lattice site 
index j. In a next step the non- interacting leads are pro- 
jected out. 12 This results in an additional diagonal and 
Matsubara frequency w-dependent one-particle potential 



0>l MeadM |jV) = Y^-V 1 + ^ 

xS 3,j' S 3,N <W 



(11) 



on site N of each wire v. The conductance of the infinite 
system Eq. (1) can then be calculated considering a finite 
system of MN + 1 lattice sites. 

The flow equations of the matrix elements with j, j ± 
1 G [1, N] are 



Ay^ A = L V V JJ A 

« 2tt ^ j ' j+r 



J= ±Ar=±l 



x (j + r,u\g {iuj) \j + r, v) 



(12) 



d _.. a u A 



^ J±1 " 2tt 



(j,v\G A (iu)\j±l,v) , (13) 



with the propagator 

g A (iu) = [g«\iu) - McadM - s A ] _1 , (w) 

which is a (MN + 1) x (MN + l)-matrix, and the initial 
condition 



^3,3 " "./,/ : I • 



(15) 
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We introduced the notation 

E# = 0>|£ A |j» • 

The matrix elements of the self-energy between the first 
sites of the wire and the dot site vanish as there is no 
interaction across these bonds. 

In a numerical solution of Eqs. (12) and (13) the flow 
starts at a large finite initial cutoff A . One has to take 
into account that, due to the slow decay of the right-hand 
side (rhs) of the flow equation for £ A , the integration 
from A = oo to A = Aq yields a contribution which 
does not vanish for Ao — > oo, but rather tends to a finite 
constant. 31 The resulting initial condition at A = A — > 
oo reads 



^3,3 

v^,A 

^3,3±1 



= (Uj-^ + U jd+1 )/2 
= . 



(16) 



As we show in the next subsection the inversion of the 
(MN + 1) x (MN + l)-matrix on the rhs of Eqs. (12) 
and (13) can be reduced to the inversion of M matrices 
of size N x N. 

At the end of the fRG flow the self-energy E A=0 
presents an approximation for the exact self-energy and 
will be denoted by S in the following. In a last step to 
obtain the Green function Q(z) which enters Eq. (8) we 
have to invert the matrix Qq 1 (z) — Vi oa d(z) — X, i.e. solve 
the single-particle scattering problem with X and the un- 
renormalizcd junction as potentials. For U ^ 0, due to 
the fRG procedure, X and thus the Green function as 
well as the conductance explicitly depend on the num- 
ber of sites in the interacting part of the wire N and the 
number of wires M. These dependences are suppressed 
in the notation we use. Similar to the case of a local- 
ized impurity in a single infinite wire 9,31,32 in each wire 
v the real space matrix elements • and £y - ±1 have an 
important spatial dependence. They show a long-range 
oscillatory behavior around an average value with an am- 
plitude which slowly decays with increasing distance from 
the junction. The scattering off this potential leads to the 
power-law scaling of the conductance. 

Considering different type of geometries of inhomogc- 
neous TLLs (single impurity, double barrier, triangular 
Y-j unction with a magnetic flux) it was earlier shown 
that the above approximation scheme leads to accurate 
results for weak to intermediate interactions such that 
1/2 < K < 3/2. 11 ' 12 ' 22 ' 26 ' 28 ' 31 - 32 In particular, in cases 
where exact expressions for scaling exponents are known 
from field theoretical models they are reproduced to lead- 
ing order in the interaction. 

Instead of analyzing the scaling of the conductance as 
a function of A for a given set of junction parameters 
as well as fixed U and N, we always integrate the flow 
equations down to A = and use the energy scale 



TTVf 



(17) 



exponents. This procedure has the advantage that each 
value of the scaling variable corresponds to a physical 
system. 



B. Scattering theory 

In this section we use single particle scattering theory 
to reach two goals, (i) We are aiming at an expression 
for the matrix elements of the A-dependent Green func- 
tion entering the rhs of Eqs. (12) and (13) that only re- 
quires the inversion of M matrices of size N x N. This 
way the numerical effort to integrate the flow equations 
can considerably be reduced, (ii) Similarly to the case 
of resonant tunneling in a single infinite wire 12 we want 
to derive equations using Eq. (8) in which the effective 
transmission is expressed in terms of the diagonal matrix 
elements of an auxiliary Green function at site 1 of each 
wire and the dot site Green function. We derive two re- 
lations of this type that will be helpful to gain a deeper 
understanding of our results. In the case of a symmetric 
junction one of them directly leads to our first result for 
the conductance. 

The Green function G(z) = [G^iz) - Vi cad (z) - X] _1 
can be understood as the resolvent matrix of an effective 
single-particle Hamiltonian h\ p (z) with a Hilbert space 
of size NM + 1. A single-particle basis is given by the 
states {\j, v) ,\d)}, where \d) denotes the Wannier state 
centered on the dot site. The single-particle version of 
the junction Hamiltonian Eq. (7) is 

M 
v=\ 

The resolvent can be decomposed as 

G(z) = Gdc(z) + Gdc{z)h iunc G(z) (19) 

with Gdc{z) being the resolvent of the disconnected wires 

G dc (z)=[z-h° lp (z)]- 1 . (20) 

The Hamiltonian h\ p follows from h\ v after taking t v = 

for all v. Applying the projector P = X)*Li Pv with 

P v = Y^j=i \ ji v )(ji v \i to tnc an d right in Eq. (19) 
one obtains 

PG(z)P = £ % + 

V 

xJ2tut^ dc (z)\l^)(l^'\Gi(z) (21) 

with 



G d (z) = 



z-^-£^<M<?dc(*)|i^) 



(22) 



as our scaling variable. It constitutes an infrared cut- 
off of any power-law scaling with interaction dependent 



and 



G'd c = PuGa c (z)P v = [zP v - P v h%{z)P v ] 1 . (23) 
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Calculating Q v Ac for fixed v requires the inversion of a 
TV x TV matrix. 

The steps leading to Eq. (21) can also be performed 
at any finite cutoff scale A. To determine the matrix el- 
ements of Q A entering the rhs of the flow equations for 
the self-energy (12) and (13) thus requires the knowledge 
of the tridiagonal parts of the inverse of M tridiagonal 
TV x TV-matrices zP v — P v h\ p {z)P u and a single column 
of each of these matrices. Numerically both can be com- 
puted in 0{N) time. 31 We can therefore easily treat a 
fairly large number M of wires each of up to TV = 10 7 
lattice sites with non-vanishing nearest-neighbor interac- 
tion. 

Along the lines of Eqs. (19) to (24) of Rcf. 12 it is 
straightforward to derive the relations 

|W(e)| 2 = 4A„(e)A^(e)|^(£ + *0)| 2 . (24) 

and 

, x,2 4A l/ (e)A y >(e) 

(25) 

with real functions £l u (e) and A„(e) given by 

Q v {e)-iA v (s)=tt(l,v\g^ c (e + iO)\l,v). (26) 

Eqs. (24) and (25) are the expressions for the effective 
transmission which can be used instead of Eq. (8). Later 
we make extensive use of these relations. 

For V = due to particle-hole symmetry f2„(0) = 
at arbitrary U. If we furthermore consider a symmetric 
junction with t\ = t<i = ... = tut = t (i^M-symmetric 
junction) all A„(0) are equal and Eq. (25) simplifies to 
\Uy(0)\ 2 = 4/M 2 independently of U and TV. The re- 
sulting conductance 

e 2 4 

with v, v' = 1, . . . , M and v ^ v' , not only follows in the 
case of a junction with t = 1, but for all t > 0. This can 
be explained as a resonance phenomenon. Within our ap- 
proximation scheme we thus obtained our first result for 
the conductance of an interacting system. As G v y is in- 
dependent of 5 n we identify the above case as a FP. More 
precisely it corresponds to a one-parameter line of FPs 
as the hopping t can be varied freely. When considering 
this case we most of the time leave the dependence on t 
implicit and refer to it as a FP. For M non-interacting 
wires 4- is the conductance maximally allowed by the 
unitarity requirement of the ^-matrix. We thus denote 
this FP as the "perfect junction" FP. As our approxima- 
tion is correct to order U we can conclude that if there is 
any interaction dependent correction to Eq. (27) it must 
be at least of order U 2 . 

To gain additional analytical insight we next study the 
one-particle spectral function for a symmetric junction of 
M wires. From the spectral function the (5,/v-dependence 
of the conductance from one of the equivalent wires to 
an additional wire that is only weakly coupled to the 
junction can be deduced. 



IV. THE HARTREE-FOCK APPROXIMATION 
FOR SYMMETRIC JUNCTIONS WITH V = 

The one-particle spectral function p bc of a TLL with 
an open boundary, taken at the chemical potential and 
close to the boundary, shows power-law scaling as a func- 
tion of 5 N with the exponent a B = l/K - 1. 4 > 6 >31,32 j^ c _ 
markably, for U > this behavior of the spectral func- 
tion can already be obtained from the non-self-consistent 
HFA, with a scaling exponent 27 



a 



HF 



U/n 
1 + U/tt ' 



(28) 



Using Eq. (6) it is straightforward to show that a§ F and 
olb agree to leading order in the interaction U. The ap- 
pearance of power-law scaling within the HFA can be 
traced back to the spatial dependence of the self-energy. 
Similar to the fRG approximation of the self-energy 31 ' 32 
the HFA self-energy shows a long-range oscillatory depen- 
dence on j that implies power-law scaling of the spectral 
weight (see below). In contrast to the boundary problem 
in the single impurity problem the HFA does not repro- 
duce the correct exponent even to leading order in U as 
the essential RG flow of the impurity is not included. 32 
One can expect that the HFA leads to meaningful results 
in all models of inhomogeneous TLLs with repulsive in- 
teraction in which such a flow is unimportant. For at- 
tractive interactions the HFA can not be used even for 
the boundary problem as oc \U\ > 0, while the exact 
exponent as is negative. We note in passing that for in- 
homogeneous TLL the application of the self-consistent 
HFA leads to unphysical results. As an artifact, the iter- 
ation of the Hartree-Fock equations generates a ground- 
state with charge-density wave order. 

The above insights motivate us to study the spectral 
function p\ on lattice site 1 (of one of the wires) for a 
symmetric junction of M wires, with t v = t for v = 
1, . . . , M and V = using the HFA. We show analytically 
that at the chemical potential 



pi oc S 



N 



(29) 

with a scaling exponent a^f that depends on M and 
U. In the following we refer to the spectral function 
evaluated at the chemical potential p — as "the spectral 
weight" . 

In the HFA and for t v = t as well as V — the Hamil- 
tonian (1) reads 

M oo 

u=l j = l 
M 

-*£( c t„ d +dci,v), (30) 



where 



t(j) 



1 + ^,j+i(cj+i,^>>° for 1 < j < TV 



otherwise 



(31) 
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with the HFA self-energy E?J +1 = U 3 . J+1 {c] +l v c J V ) . 
The expectation value (. . .) is taken with the many-body 
groundstate of H kin + H ]unc (V = 0) = Hn F (U = 0). As 
we shifted the density by its average value the Hartree 
term vanishes. 

The normalized single-particle eigenstates \^ e , I) of the 
one-particle version /ihf of the Hartree- Fock Hamiltonian 
can be classified according to their behavior in 2ir/M- 
rotations (^M-symmetric case) which leads to the ex- 
pansion 

M oc 

|* £ ,0=^oa ( £ )M) + ^^ e< ^ I/ 4°( e )lj» ' ( 32 ) 

v=lj=l 

with coefficients a^p(e) and a (e). For odd values of M 
one has I = 0, ±1, . . . , ±(M — l)/2, while for even values 
Z = 0,±l,...,±M/2-l,M/2. 



A. The eigenstates for [7 = 

In a first step we determine the eigenstates of the 
non-interacting system. The energies are e = e(k) = 
— 2cos(fc). For fixed energy the wave number is thus 
given by k(e) — arccos(— e/2). Using Eq. (32) 



/iHF(tf = 0)|tf e) Z>o =e|* e ,0o 



yields 



The resulting single-particle states \^ £ , Z)o can then be 
used to calculate the groundstate expectation value that 
enters the HFA Hamiltonian 



( c j+l,i/ c j,v)o 



E 



deaf +1 {e)a>;>{e) 



(0/ 



B/2 



—Y 



dk {cos[fc] 
cos [k(2j + 1) + 2Si(k)}} 



(39) 



with the Fermi wave number kp — ir/2 and band width 
B = 4. As described below to determine the scaling of 
the spectral weight of the interacting system for <5jv — * 
we only need to know the behavior of this expectation 
value for 1 <C j < N 



v _ 1 



1 - 



til. 

2j + l 



1 -^) + °(^ 2 ) 



This asymptotic behavior can be obtained from Eq. (39) 
using integration by parts and S (kp) = ±ir/2. For I < 
j < N — j s the interaction Ujj+i Eq. (4) takes its bulk 
value U and it follows that 



7T 

- U + 2 

7T 



7T 2j + 1 



2 



1 ^HF 

I - a R 



-iy 



2j + 1 



,(40) 



for 1 < j < AT - j s 



-tSi t0 Sj t iao(e) (33) 

for the coefficients of |* e , Z) . For Z = there is an addi- 
tional equation 

ea (e) = -Mta { °\e). (34) 
The solution of Eqs. (33) and (34) is 

af{e) = A(M, e) sin {j k{e) + Si [k(e)}} (35) 

and 

ao(e) = A(M, e) sin {<5 [k(e)]} /i , (36) 
with the phase shifts 

j iW = [arctan{^tan(fc)} for Z = (3?) 
[ otherwise , 



the normalization factor 



A(M, e) = 



nMv(e) 

and the velocity v(e) =2 sin [k (e)]|. 



(38) 



B. Spectral weight for U ^ 

Within the HFA the spectral weight on the first site of 
one of the equivalent legs is determined by the amplitudes 
of the e = eigenstates |\& e =o, Z) of /ihf 



Pl cx^|(I,z,|vI/ £=0 ,Z)| 2 =^ a«(0) 



(41) 



To avoid proliferation of symbols or indices the expan- 
sion coefficients of the interacting HFA eigenstates are de- 
noted by the same symbols af 1 and ciq as the coefficients 
of the non-interacting eigenstates in the last subsection. 
With the expansion Eq. (32) the stationary Schrodingcr 
equation 



/ihf|* £ =o,Z) = 



leads to coupled equations for the coefficients efp (0). The 
Z = equation (34) also holds for U 7^ which for e = 

leads to a^ 0) (0) = 0. Only the I = eigenstate |* e =o,0) 
has a non- vanishing amplitude on the dot site which im- 
plies that the af\o) for Z ^ can be calculated as for a 
semi-infinite chain. For j > 2 and s = the Schrodingcr 
equation gives 

(3, ^hf|* £ =o,0 = = -t(j>?+i(0) - t(j - 1)4^(0) . 
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This relation can be solved iteratively leading to 



4? + ,=(-D j 4"(o)n^ 



t(2i -1) 
0~ 



(42) 



Because of 

(l,l/|/»HF|* e =0,i> = 0=-t(l)4°(0) , 

a^(0) = for all even j. 

Without loss of generality we now consider the case of 
even TV. As the interaction Ujj+i vanishes for j > N, 



t(2i-\) = t(2i) for i > N/2+1 and 



hj+i 



is independent 



of j for j > N/2 + 1. Together with the asymptotic 

scattering state form of the a,j (0) Eq. (35) this implies 
that with respect to the ^-dependence we find 

2 



pi oc yjq 



oc 



N/2 

n 



t(2i) 



L = l t(2i - 1) 



(43) 



first lattice site pi and the spectral weight on the dot 
site p<i are inversely proportional to each other. It thus 
follows that 



Pd oc 5 



N 



(46) 



For M = 1 the dot site is the last site of a semi-infinite 
chain with open boundary conditions and a^ F agrees 
with a§ F . For M = 2 the dot site corresponds to a bulk 
lattice site of a homogeneous TLL. As long as V = 
this holds even for i ^ 1 due to a resonance. In a ho- 
mogeneous TLL the spectral weight scales as S N , with 
a = (K^ 1 + K — 2) / 2} Replacing K in this expres- 
sion by the leading order term Eq. (6) it follows that 
a = 0(U 2 ). The HFA exponent can only expected to 
be correct to leading order in U . Consistently we find 
0. At M = 2 and U > the sign of a^f changes 



a 



HF 



from aff > for M = 1 to o& F < for M > 3. 



M 



The conductance across a weak link 



We next evaluate the rhs of this expression for large N. 
We do this separately for the numerator and denomina- 
tor. Using Eq. (40) one obtains 



In 

N 



N/2 

IJt(2i) 

i=i 



N/2 

E 

i=l 



In [t(2i)} 



■ a 



HF 



1 



M 



\nN + O{N ) . 



The first two terms follow from the factors of the product 
in which t(2i) can be replaced by Eq. (40). The remaining 
factors lead to the last summand of order N°. Similarly 
one gets 



In 



N/2 

IJt(2i-l) 



i=l 



N/2 

= £>[t(2i-l)] 



i=l 



= — In H 

2 V 7T 



l-jj)lnN + O(N ). 



Combining these two asymptotic expressions leads to the 
power-law scaling of the spectral weight on the first lat- 
tice site of each wire 



pi oc N" M oc S 



N 



with 



^hf 



^HF / ^ 



U/tt 
1 + U/n 



(44) 



(45) 



This derivation shows that the power-law scaling of the 
spectral weight directly follows from the long-range spa- 
tial oscillations of the self-energy. 

Using Eqs. (21) and (22) it is straightforward to show 
that for a symmetric junction the spectral weight on the 



We next consider a junction of M wires, M — 1 of 
them with hopping t from the first lattice sites of the 
wires to the dot site, while the leg with v — M is cou- 
pled by the hopping amplitude tu *C t. The HFA result 
for the spectral weight can be used to determine the scal- 
ing exponent of the conductance from one of the M — 1 
equivalent wires to the wire M. Applying Fermi's golden 
rule one can argue that in this tunneling limit the scaling 
of the conductance across the weak link is determined by 
the product of p Q b c and pd leading to 



G v ,m oc S 



N 



for v 



1, ... ,M - 1, with 

7 HM) = 2ag F /(M - 1) . 



(47) 



(48) 



This constitutes our second result for the transport 
through a dot junction. For M = 2 the junction prob- 
lem is equivalent to the single impurity problem in the 
limit of a weak link (strong impurity) which is charac- 
terized by the exponent 2as. 4 We reproduce this result 
to leading order in U. Anticipating the RG language 
Eq. (47) indicates that for U > a weak link of a sym- 
metric (M — l)-leg junction to an additional wire is an 
irrelevant perturbation. The exponent 7^ F (M) > 0, i.e. 
G v ,m — > and with decreasing <5/v the system "flows 
back" to the "perfect junction" FP of the (M - l)-wire 
system with the conductance Eq. (27). 

Without using the golden rule like argument and ap- 
plying the fRG we next numerically confirm the power- 
law scaling of G Vi m as well as the RG interpretation of 
the HFA results. Going beyond the symmetric junction 
with V = and an additional weak link, we study the 
conductance for general junction parameters. The fRG 
can also be used for U < 0. 
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V. FRG RESULTS 

In this section we present the results for the FPs and 
the scaling of the conductance obtained from numerically 
solving the fRG flow equations of Sect. Ill for U ^ 0. 
In subsections VA and VB we investigate two specific 
classes of junction parameters. The results from these 
cases can be combined and lead to the comprehensive 
picture for arbitrary junction parameters presented in 
subsection VC. 



A. A symmetric junction with one modified link 

We here consider a junction of M > 2 wires, M — 1 
with hopping i (not necessarily t = 1) between the first 
site of the wire and the dot, while the additional one has 
hopping tM =/= t. The dot site energy V we set to zero 
and the real part f2„(0) of the auxiliary Green function 
Eq. (26) vanishes. In addition A„(0) = A„,(0) = A for 
all u, v' < M — 1. Two cases can be distinguished. 



1. A weak link 

The first one is the weak link situation already treated 
by the HFA with r = tu/t -C 1. We thus slightly perturb 
the "perfect junction" FP of a (M — l)-wire system in a 
specific way. As expected we numerically find that the 
effective transmission between the M — 1 legs with t is 
close to the perfect transmission 4/(M — l) 2 . For U > 
and decreasing Sn the fRG data for \t v y | 2 approach this 
value, while they leave it for U < 0. At U = 0, \t v y\ 2 
is independent of 8n- Similarly |i„,M| 2 — * for U > 0, 
while it takes a small increasing value for U < 0. We 
thus analyze the power-law scaling of |£i/x| 2 and |^,m| 2 
as a function of Sn with respect to 4/(M — l) 2 and 0, 
respectively. The effective exponents as a function of 5n 
obtained by taking the log-derivative of the fRG data 
calculated for M = 5, t = 1, tu = 1CP 3 , and different U 
is shown in Fig. 2 on a log-linear scale. For small 5n and 
all U both exponents approach the same [/-dependent 
plateau value, which is our fRG approximation for the 
scaling exponent. Even for larger arguments the two Sjy- 
dependent effective exponents are indistinguishable on 
the scale of the plot. 

The M-dependence of the scaling exponent for differ- 
ent U is shown in Fig. 3 on a reciprocal-linear scale (sym- 
bols). The data can be fitted by (lines) 

7i(M)=/3 s /(M-l), (49) 

with f3 s being the [/-dependent fRG approximation for 
2a b = 2(1 /K — 1) obtained in Rcf. 12 for a weak link 
(strong impurity) in an infinite wire, i.e. the M — 2 case. 
The exponent (3 S agrees with 2ots to leading order in U. 
It has higher order corrections which bring it close to 
2cs.b even for intermediate U. For the four interactions 




FIG. 2: Effective exponents of pjzip ~ \^v.v'\ Z (circles) and 
\t v ,M\ 2 (stars) as a function of 8n for M = 5, t = 1, tu = 
10 , and different values of U = — 1, —0.5, 0.5, 1 from top to 
bottom. On the scale of the plot the results obtained from 
ngzrp — \t v y | 2 and |i„,M| 2 are indistinguishable. The scaling 
exponent is read-off when a plateau is reached for Sm < 5 • 
10~ 4 , which roughly corresponds to N > 10 4 . 




1 2 3 4 8 100 

M-1 



FIG. 3: Scaling exponent of the conductance for a weak link as 
a function of M — 1 (symbols) for different U. Note the recip- 
rocal scale of the x-axis. The lines show 71 (M) = (3 S /(M — 1). 



U = —1, —0.5, 0.5, 1 the value for j3 s can be read off from 
Fig. 3. At U = 1 the relative difference between the exact 
exponent obtained from Eq. (6) and our approximation 
is roughly 5%. A detailed comparison of j3 s and 2as is 
given in Fig. 5 of Ref. 12. To leading order in U Eq. (49) 
confirms the result deduced from the combined use of the 
HFA and Fermi's golden rule arguments Eq. (48). For 
M — 2, the fRG approximation (3 S is generically closer 
to the exact exponent 2a b than 2a§ F . 

For U > 0, 71 > and the "perfect junction" FP of 
the (M — l)-wire system is stable towards weakly cou- 
pling an additional wire, while for U < 0, 71 < and 
the FP is unstable. In the latter case the system ef- 
fectively incorporates the weakly coupled leg and flows 
to the "perfect junction" FP for M wires. For small to 
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>4/16 




10 10 



FIG. 4: One-parameter scaling plot of the effective transmis- 
sions |i„,M| 2 (lower curve) and \t„y\ 2 with v,v' < M — 1 
(upper curve) for M = 4, U = —1, and t = 1. The variable 
is Sn /So(tM), with a non-universal scale 5o(tAf). The asymp- 
totic values of the transmission 4/(M — l) 2 and reached for 
Sn/So^m) — ► oo as well as 4/M 2 reached for 6n /#o(£m) — ► 
are indicated on the j/-axis. 



intermediate \U\ and starting close to the "perfect junc- 
tion" FP of the (M — l)-wire system, i.e. for tu <C t, 
exponentially small 8n are required to reach the new 
FP. Even though we can treat very large TV such small 
scales are beyond the possibilities of our method. Sim- 
ilar to the single impurity problem 4, 10-12 the flow from 
one to the other FP can be shown considering data sets 
(Sn, \t».u'\ 2 ) obtained for different tu €]0, 1[ using a one- 
parameter scaling ansatz. For fixed U the conductance 
is only a function of Sn/Sq^m), with an appropriately 
chosen non- universal energy scale <5o(£m) and the differ- 
ent data sets can be collapsed on a single curve. This is 
shown in Fig. 4 for the two different effective transmis- 
sions I^mI 2 (lower curve) and \t v y | 2 with v, v' < M — 1 
(upper curve) and M = 4, U = — 1, f = 1 on a log- 
linear scale. It is important to note that the effective 
transmission 4/M 2 between all legs is not achieved by 
an increase of tu- As there is no interaction across the 
links of the dot to the first sites of the wires, the hopping 
across these links is independent of the RG cutoff (here 
Sn)- The transmission 4/M 2 follows from the build up 
of a self-energy during the RG flow, that, interpreted as 
an effective scattering potential, leads to a resonance at 
the chemical potential. For 71 < 0, that is U < 0, the 
plateau reached in Fig. 2 does not present the asymptotic 
behavior for Sn — > and considering even smaller 5n a 
deviation from the plateau value 71 can be observed. 

We next show that the scaling behavior of the conduc- 
tance can be understood from the scaling of the spectral 
weights on the first site of one of the equivalent legs and 
the first site of the additional wire. For e = we get 
from Eq. (25) 



\tv,v' | 2 = 



(M - l) 2 (M - l) 3 



D + 0(t 4 ) 



(50) 



for v ^ v' and v, v' < M — 1 , as well as 
4 



\tv,M\ 



D z + 0(r b ), (51) 



(M-l) 2 (M-l) 3 " 
for v < M - 1, with 

D = A M (0)/Aocr 2 . (52) 
For the reflection back into wire v < M — 1 it follows 

Rv = 1 — ^ \~tv,v'\ 2 

(M - 3) 2 4 



+ 



(M-l) 2 (M-l) 2 
D 2 + 0{t 6 ) . 



(M-l) 



D 



(M- l) 3 



(53) 



As discussed in Sect. Ill the auxiliary Green function is 
calculated with a self-energy which has been determined 
in the presence of the weak link. Thus A/i 2 and the 
spectral weight p\ of a perfect (M — l)-leg junction differ 
to order r 2 . If corrections of order r 2 are neglected in A 
and Am(0), D can be replaced by p bc/pi- We numer- 
ically find that for M = 1,2,... the spectral weight on 
the first site of a perfect M-leg junction scales as 



Pi cx S 



with 



N 



M 2 



2 

M 



(54) 



(55) 



consistent with the derivation using the HFA in Eq. (45). 
The above replacement thus leads to D cx <5^ < ' M ~ 1 \ 

3 /2 

where we used p Q b c oc 8^' . Inserting this expression in 
Eqs. (50) and (51) we reproduce the results directly ob- 
tained from calculating the effective transmission. Note 
that for M — 3 the second term in Eq. (53) cancels and 
the reflection scales with an exponent that is twice as 
large as the one of the transmission. Further down we en- 
counter another case in which the prefactor of the leading 
order term vanishes for specific parameters, which leads 
to a doubling of the scaling exponent. 

In the limit M — > 00 the dot site is coupled to so 
many legs that it creates an infinite barrier. Consistently 
Eq. (55) gives limM^oo a 



fRG 

M 



-4r and the spectral 



weight on the first sites scales as the weight next to an 
open boundary. 



2. A slightly modified link 

As the second example with a single modified link we 
consider |*m — i|/£ = t <C 1. We now analyze the scaling 
of the effective transmission within the subsystem of the 
M — 1 equivalent wires and into the leg with the modified 
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3 3 4 8 101 

M 



FIG. 5: Scaling exponent of the conductance for a slightly 
modified link as a function of M for different U. Note the 
reciprocal scale of the x-axis. The lines show the fit Eq. (56). 



hopping with respect to the transmission 4/M 2 of the 
perfect case. Similarly to Fig. 2 exponents can reliably 
be extracted for sufficiently large N. 

For M > 3 the exponents of |i„.„<| 2 and I^m) 2 , with 
v, v' 6 [1, M — 1] turn out to be equal. For several U the 
M-dependence of the scaling exponent is shown in Fig. 5 
on a reciprocal-linear scale. To first order in U it is given 
by —(3 S /M. More accurately the data can be fitted by 



and (51) (which is valid for all M > 2) 

2 4 4(J>-1) / 2\ 

with D as defined in Eq. (52) . The difference D — 1 shows 
power-law scaling with exponent 72 (M). For M = 2 
the prefactor of the leading order term linear in D — 1 
vanishes and the scaling exponent of I^mI is doubled. 
Inserting M = 2 in Eq. (56) gives 272 (2) = (3 W as the ex- 
ponent characterizing the deviation from transmission 1. 
This result is expected since the case M — 2 corresponds 
to the situation of a perfect infinite wire interrupted by 
a small impurity characterized by the scaling exponent 

In the single impurity problem the fRG approximations 
fl s for 2a B = 2(1 /K — 1) and j3 w for 2(K — 1) are correct 
to order U. To leading order 2(K - 1) w -2(1/ K - 1) 
[see Eq. (6)] and the second term in Eq. (56) is of order 
U 2 . As terms of order U 2 are only partly included in our 
approximation scheme it is questionable if a term similar 
to the second one will be present in the (unknown) exact 
expression for 72 . To derive the dependence of 71 and 72 
on K and M presents a challenge for any method which 
does not require approximations in the strength of the 
interaction. 



72(M) = ~M + 2 ^r- 



(56) 



with the exponent (3 W found in Ref. 12 for the scaling 
of the transmission through a local weak impurity. To 
leading order in U, /3 W agrees with the exact weak single 
impurity exponent 2(K — 1). A comparison of 2(K — 1) 
and our fRG approximation is given in Fig. 7 of Rcf. 12. 

For U > 0, 72 < and the "perfect junction" FP of the 
M-wire system is unstable towards changing one of the 
hopping amplitudes. For tM < t the system flows to the 
"perfect junction" FP of M — 1 wires. This follows from 
the build up of a self-energy during the RG flow, that 
leads to a vanishing transmission from wire v < M — 1 
to wire M, but to a resonance with perfect transmission 
between the M — l equivalent legs. For tM > i the flow is 
to a FP with vanishing conductance between all wires - 
the "decoupled chain" FP. The vanishing of the conduc- 
tance follows from the long-range oscillatory behavior of 
the self-energy generated in the RG flow. In both cases 
the flow from one to the other FP can again be shown 
using a one-parameter scaling ansatz. For U < 0, 72 > 
and the "perfect junction" FP of the M-wire system is 
stable, regardless of the sign of t — tM- 

For M = 2 there is no \t„y \ 2 and 1 — \t u ^ M | 2 scales with 
272(2). The appearance of the factor 2 can be explained 
by considering an expansion of \t u .M\ 2 similar to Eqs. (50) 



B. A symmetric junction with a dot site energy 

As our second specific case we study a symmetric M- 
leg junction with t v = i for v = 1, . . . , M and a non- 
vanishing dot site energy V . Then, due to symmetry, all 
fi„(0) and Aj,(0) [see Eq. (26)] are equal and we suppress 
the index v. Eq. (25) simplifies to 



l*| 2 = 772 



M 2 A 2 (0) 



M 2 [V + Mft(0)] 2 + M 2 A 2 (0) 



(58) 



The transmission is determined by the single complex 
parameter 



g = V + Mfi(0) - iMA(0) = 



Gd(0 + i0) 



(59) 



which itself is a function of the junction parameters (t, V) 
and M. Via the RG flow of the self-energy for U ^ it 
moreover develops a dependence on Sn and the interac- 
tion U. The RG flow can nicely be visualized by plot- 
ting g in the complex plane with Sn as a parameter. 12 ' 22 
This is done in Figs. 6a) to 6c) for M = 2, 3, 10, V = 
0.5,1,1.5,2, and different t e [0.01,0.7]. For decreasing 
8 N e [5 • 10~ 4 , 2.5 • 10- 1 ] and U ^ each fixed param- 
eter set (t, V) leads to a flow line. The general form of 
the flow diagrams is independent of the absolute value of 
U. The data of the figures were calculated for U = —1, 



11 




FIG. 6: Flow of g Eq. (59) for the case of a symmetric junction with on-site energy V = 0.5, 1, 1.5, 2, different i € [0.01, 0.7], 
M = 2, 3, 10 legs, and U = —1. The arrows indicate the direction of the flow for U < 0. For U > it is reversed. 




-0.5 



4 
M 



• U=-1 

* U=-0.5 
■ U=0.5 

♦ U=1 



100 



self-energy, which leads to a a resonance at the chemical 
potential. For U > the flow is clockwise and the line of 
"perfect junction" FPs is unstable. The system flows to 
the line of stable "decoupled chain" FPs. The vanishing 
of the conductance is a consequence of the long-range os- 
cillatory dependence of the self-energy on j. For M > 3, 
U < 0, and 6n — * 0, -hag diverges while Reg goes to a 
constant, which implies that the system reaches a "per- 
fect junction" FP. For U > all trajectories approach the 
x-axis at V and the conductance vanishes (line of "de- 
coupled chain" FPs). To determine scaling exponents we 
next separately consider small and large V. 



FIG. 7: Scaling exponent of the conductance for a small on- 
site energy as a function of M (symbols) for different U. 
Note the reciprocal scale of the a;-axis. The lines show the 
fit Eq. (60). 



which leads to the flow direction indicated by the arrows. 
For U > the direction is reversed. On the line Im g = 
the conductance vanishes and the cc-axis forms a line of 
"decoupled chain" FPs. For Keg = the transmission is 
4/M 2 and all points on the y-axis are "perfect junction" 
FPs (line of FPs). The "perfect junction" transmission 
4/M 2 is also reached at g — oo provided that this point 
is approached such that — Img — > oo and Reg goes to a 
constant. 

For M = 2 the flow approximately follows a section of 
a circle centered around the origin with a radius V. For 
U < the direction is counter-clockwise and the line of 
"perfect junction" FPs (y-axis) is stable. The on-site en- 
ergy V does not get renormalized in the RG flow as the 
interaction between the dot site and the first sites of the 
wires is assumed to be 0. The perfect transmission 4/M 2 
does thus not follow from a decrease of V during the RG 
flow but is a consequence of the spatial structure of the 



1. A small on-site energy 

For small V, the transmission between all the equiva- 
lent wires is close to the perfect value 4/M 2 . The depen- 
dence on Sn can be described by a power-law. In Fig. 7 
the scaling exponent as a function of M is shown for dif- 
ferent U. To leading order in U it is independent of M. 
The exponent can be fitted by 



73 (M) = -/? s + 4 



flu 



M 



fls ( _ J_ 

M 



(60) 



In accordance with the stability properties of the line 
of "perfect junction" FPs discussed in connection with 
Figs. 6a) to 6c), 73 < for U > and 73 > for U < 0. 
The M = 2 case is equivalent to the problem of a single 
weak site impurity interrupting an otherwise perfect in- 
finite chain and 73(2) is equal to the respective scaling 
exponent fl w (similar as for a slightly modified link in 
Sect. VA). Within our approximation scheme the sec- 
ond term in Eq. (60) is thus important to reproduce a 
result obtained earlier for M = 2. As it is of order U 2 
it is nonetheless unclear if a similar term occurs in the 
(unknown) exact expression for 73. 
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2. A large on-site energy 

In the limit of large V the transmission between all the 
wires is small and we analyze the scaling with respect to 
zero transmission. The scaling exponent is independent 
of M and up to our numerical accuracy given by 



74 = (3 S 



(61) 



i.e. the strong impurity (weak link) exponent for a single 
infinite wire (see Fig. 5 of Ref. 12). 



C. Fixed points and renormalization group flow for 
general junction parameters 

The RG flow, FP structure, and scaling exponents for 
arbitrary junction parameters can be understood based 
on the results obtained for the above two classes of junc- 
tion parameters. To shorten the discussion we here focus 
on the more important case U > 0. Results for U < 
can be deduced by inverting the direction of the flow. 

As found in subsections V A and V B the "perfect junc- 
tion" FP of the M-leg system with effective transmission 
4/M 2 is unstable towards the two possible perturbations 
in which only one junction parameter is modified: (i) a 
single modified hopping and (ii) a non-vanishing on-site 
energy. These instabilities are characterized by the two 
exponents 72 (M) (modified link) and 73 (M) (on-site en- 
ergy). 

In the case (i), for M > 3, and t larger than the sin- 
gle modified hopping the system flows to the "perfect 
junction" FP for M — 1 legs. For Sn — > the conduc- 
tance across the modified link vanishes with exponent 
71 (M), while the conductance within the subsystem of 
the M — 1 equivalent wires approaches 4/(M — l) 2 with 
the same exponent. For M = 2 the flow is to the "de- 
coupled chain" FP and the conductance vanishes with 
exponent 71(2) = 74 = f3 s . If the hopping between the 
first site of wire v = M and the dot site is larger than t 
for S N and all M > 2 the "decoupled chain" FP is 
reached. To analyze the scaling of the conductance close 
to this FP we start out from Eq. (25) and use an expan- 
sion similar to Eqs. (50) and (51). The small parameter 
of the expansion is D = A/ Am- The numerics shows 
that it asymptotically vanishes as 6^ . For v < M — 1 
this leads to 



t u are reduced compared to the hopping t across the re- 
maining M2 = M — Mi links the system flows to the 
"perfect junction" FP of the M 2 -leg system. If Mi of 
the t v are increased compared to t the system generically 
approaches the "decoupled chain" FP. This behavior is 
changed if M\ > 2 of the M\ increased hoppings are equal 
and larger than all the others. In this case the system 
flows to the "perfect junction" FP for M\ legs. If all t v 
are different the "decoupled chain" FP is reached. 

If the perturbation of the "perfect junction" FP of the 
M-leg system consists of an on-site energy V and one or 
more modified hoppings the system flows to the "decou- 
pled chain" FP. 

From these considerations we conclude that for the 
majority of possible junction parameters of the M-wire 
junction, for 8m — ► the system flows to the "decoupled 
chain" FP. On asymptotically small scales the conduc- 
tance vanishes as 



G v v i oc S^? 3 . 



(64) 



Depending on the junction parameters as well as on the 
wire indices v and z/, £ might be 1 or 2 (for examples, 
see above). 

Only for V = and if M max > 2 of the M links between 
the dot site and the first lattice sites of the wires have 
hopping i max , where i max = max„ e [ ljM ]{i„}, the system 
flows to the "perfect junction" FP for M max legs. In this 
case the power-law scaling of the conductance between 
two wires coupled to the dot with hopping i max is given 

by 



h ML 



G v y cx 8 N 



Ps/M„ 



(65) 



The conductance between one wire coupled by i max and 
the other by t v < t max vanishes as 



G v y CtO N ' 



while it goes like 



G v y oc 5 2 ^ S ^ M " 



(66) 



(67) 



if both wires v and v' have a hopping to the dot site 
smaller than i max . As in Eq. (63) the factor 2 in Eq. (67) 
follows from an expansion similar to the one used in 
Eqs. (50) and (51). 



G vM oc (62) 

while for v, v' G [1, M — 1] we find 

G vy oc 4 74 . (63) 

In the case (ii) the RG flow also ends at the "decoupled 
chain" FP and all conductances vanish with exponent 74. 

We next perturb the "perfect junction" FP of the M- 
leg system by slightly modifying more than one of the 
hopping matrix elements. If 2 < M\ < M — 1 of the 



VI. SUMMARY AND PERSPECTIVES 

In this paper we studied the conductance G v y for a 
dot junction of M semi-infinite quantum wires. The junc- 
tion as well as the wires are described by a microscopic 
lattice model. In a finite section of N sites the fermions in 
each wire are modeled to interact via a nearest-neighbor 
interaction which smoothly vanishes close to the lattice 
sites j — N (smooth contact). Investigating the scal- 
ing of the conductance as a function of Sn oc 1/N led 



13 



to a comprehensive picture of the FP structure, scal- 
ing exponents, and RG flow of the model studied. We 
mainly used an approximation scheme that is based on 
the fRG method and provides reliable results for week to 
intermediate interactions with 1/2 < K < 3/2. At half- 
filling this corresponds to bulk nearest-neighbor interac- 
tions — 1.5 < U < 2. Additional insights were obtained 
using the HFA. 

Compared to the well studied single impurity problem 
for M > 3 the low-energy physics of the M-wire junc- 
tion is much richer, allowing for a variety of FPs and 
scaling exponents. Furthermore for M > 3 one of them, 
the "perfect junction" FP, is characterized by two ex- 
ponents: 72(M) and ^(M). They can individually be 
read off from the scaling of the conductance if the FP 
is perturbed in a specific way. In contrast to the sin- 
gle impurity problem, in which the FP reached is solely 
determined by the sign of the interaction, for junctions 
of three and more wires this in addition depends on the 
junction parameters. Depending on the wire indices be- 
tween which the conductance is calculated the scaling ex- 
ponents close to a FP might differ by factors of two. This 
can have two reasons: (i) In an expansion of the effective 
transmission in terms of a small parameter that carries 
the power-law scaling [see e.g. Eqs. (50) and (51)], de- 
pending on v and v' the first or second order term might 
be the first non-vanishing one. (ii) For certain M the 
prefactor of the leading order term might vanish. 

Junctions and networks of TLLs were earlier 
studied. 15-21 The method and model used in Ref. 17 
come closest to ours. The authors apply a fermionic 
poor-man's like RG originally developed for the single 
impurity problem to the three- and four-leg junction. 9 
They consider a microscopic dot junction model to mo- 
tivate the investigation but then leave the framework of 
this model when studying the RG flow of an effective 
S'-matrix. The results obtained in this paper for three 
legs are partly equivalent to ours. In Refs. 12 and 26 a 
detailed account of the differences of the poor-mans RG 
and our method is given for the single impurity case and 
resonant tunneling in a TLL. 

Using the fRG and HFA we obtained expressions for 
the scaling exponents which we expect to be correct (at 
least) to leading order in U. It is very desirable to derive 
the exact K- and M-dcpcndence of these exponents. For 
increasing complexity of the junctions obtaining such ex- 



pressions for simplified, effective, but still generic models 
requires very sophisticated methods, even if only specific 
parts of the parameter space are considered. 19 

We considered a dot junction model in which the in- 
teraction between fermions on the dot site and the first 
lattice sites of the wires is set to zero. We also investi- 
gated the case in which this interaction does not vanish. 
The additional nearest-neighbor interaction across the M 
bonds does not alter the results presented here. 

An extension of our method to the case in which the 
strength of the bulk interaction depends on the wire in- 
dex v is straightforward and might lead to new insights. 
As exemplified in Ref. 22 the fRG can also be used to 
investigate other types of junctions (e.g. ring like geome- 
tries), with different FPs and new scaling exponents. 

For the complex junction studied here the temperature 
T and the infrared cutoff Sn present equivalent scaling 
variables only on asymptotically small scales. The fRG 
method can be set up for finite T and leads to reliable re- 
sults also for intermediate to large temperatures. 12 ' 26 For 
two wires our model is equivalent to the one considered to 
study resonant tunneling through a quantum dot embed- 
ded in a TLL, with a one lattice site dot. In this case, for 
fixed N, and 5n < T < B the conductance as a function 
of T shows a very rich behavior. Depending on the dot 
parameters, temperature regimes in which G(T) follows 
"universal" power-laws as well as non-universal regimes 
were identified. 12,26 We thus expect to find similar rich 
behavior for M > 3. The conductance as a function of 
temperature is easily accessible in transport experiments. 
Investigating G v y (T) might thus also become important 
for the interpretation of future transport experiments on 
junctions of quasi one-dimensional quantum wires. We 
will present results for G v y (T) in an upcoming publica- 
tion. 
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